Confinement-induced Berry phase and helicity-dependent photocurrents 
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The photocurrent in an optically active metal is known to contain a component that switches 
sign with the helicity of the incident radiation. At low frequencies, this current depends on the 
orbital Berry phase of the Bloch electrons via the "anomalous velocity" of Karplus and Luttinger. 
We consider quantum wells in which the parent material, such as GaAs, is not optically active 
and the relevant Berry phase only arises as a result of quantum confinement. Using an envelope 
approximation that is supported by numerical tight-binding results, it is shown that the Berry phase 
contribution is determined for realistic wells by a cubic Berry phase intrinsic to the bulk material, 
the well width, and the well direction. These results for the magnitude of the Berry-phase effect 
suggest that it may already have been observed in quantum well experiments. 
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The "Fermi liquid" theory of metals developed by Lan- 
dau captures the microscopic details of band structure 
and interactions in a small number of parameters that 
quantify how the quasiparticles in the metal differ from 
free electrons. While this theory has been successful in a 
wide variety of materials, it is now considered to be in- 
complete: it does not include physics resulting from the 
Berry phase of the Bloch wavefunctions. Berry phases 
result when a quantum-mechanical wavefunction changes 
smoothly as a function of some parameter pQ . The orbital 
Berry phases in metals cause an "anomalous velocity" of 
Bloch electrons, derived by Karplus and Luttinger [2] and 
later interpreted as a Berry-phase effect [3j H] . 

The Berry phase that influences basic transport prop- 
erties vanishes in materials that are symmetric with re- 
spect to both inversion and time-reversal. In metals that 
break time-reversal (TR) symmetry (ferromagnets and 
antiferromagnets) , the Berry phase leads to an "intrin- 
sic" mechanism of the anomalous Hall effect (AHE). The 
other broad class of materials with nonzero Berry phase 
are those in which TR is preserved but inversion symme- 
try is broken; in insulators, these Berry phases underlie 
the modern theory of polarization [5j [6]. Berry phases 
from spin-orbit coupling lead to "topological insulators" 
in two [9] and three [10] HQ E] dimensions. 

By contrast, at present there is no experimental obser- 
vation that has been associated with the non-vanishing 
Berry phase that is expected in metals that break in- 
version symmetry. However, the helicity-dependent pho- 
tocurrent in an optically active metal contains a Berry- 
phase contribution, as shown recently by Deyo et al. [15] , 
who give a semiclassical transport analysis of photocur- 
rents at linear order in applied intensity. We believe that 
this effect, also known as the circular photo-galvanic ef- 
fect (CPGE), is quite fundamental and differs in impor- 
tant ways from previous Berry-phase phenomena: it is 
nonlinear, frequency-dependent, and controlled by the 
Berry-phase contribution at low frequency. 

In this paper we compute the Berry-phase contri- 



bution to helicity-dependent photocurrents in realistic 
circumstances, in order to allow quantitative compari- 
son with recent experiments on semiconductor quantum 
wells. In these systems, the Berry phase contribution 
vanishes in the bulk; the the photocurrent is generated 
by quantum confinement. We first give a simple deriva- 
tion of the Berry-phase contribution in two dimensions in 
the relaxation-time approximation. A microscopic tight- 
binding calculation for a model GaAs (110) quantum well 
is then described. It yields a magnitude of the photocur- 
rent that agrees well with experiments by the Regensburg 
group p~4j [15] [16]. Finally, we show via an envelope- 
function approach that the confinement-induced Berry 
phase in a general zincblende quantum well can be pa- 
rameterized in terms of only two numbers, one intrinsic 
to the bulk material and the other a geometrical factor 
determined by the surface orientation and well width. 

Our starting point is the semiclassical equation for a 
wavepacket of Bloch electrons. The velocity vector has 
two terms [2] [5] [17] . One is the familiar semiclassical ve- 
locity, which is determined near a parabolic minimum by 
the effective mass, and the other results from the change 
in the spatial location of the electron within the unit cell 
as its wavevector k moves through the Brillouin zone: 



dx a = 1 de n (k) T ab (]) dk b 
dt ft dk a n [ } dt' 



(1) 



In the second term, known as the anomalous velocity, the 
Berry flux is written with the same symbol T ah as the 
electromagnetic field tensor in order to stress certain sim- 
ilarities. (The band index n is suppressed from now on 
as we assume a single spin-degenerate partially occupied 
band). T ah is the curl of a "vector potential" obtained 
by momentum-space derivatives of the periodic part of 
the electronic eigenstates Mfc(x) (here d a = dk a )'- 



f ab (k) = d a A b (k)-d b A a (k), 
A a (k) = -i(u k (x)\d a u k (x)). 



(2) 
(3) 



As in electromagnetism, A a is gauge-dependent (i.e., 
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FIG. 1: Snapshot of electron distribution in a circularly po- 
larized wave. The electron distribution is pushed away from 
the ground state (solid circle) slightly for small cot. At the 
time shown, the conventional velocity vo and anomalous ve- 
locity vi are in opposite directions. After a half-period, the 
conventional velocity changes sign but the anomalous velocity 
does not. The directions assume ft — f3k x z with (3 > 0. 



changes under a phase change of Bloch states) but T ah is 
gauge-independent. It is also antisymmetric, and time- 
reversal symmetry requires ^ a6 (k) = — T ah {— k). 

In two dimensions, the Berry flux has only one nonzero 
component T 12 = — JT 21 , which can be pictured as a vec- 
tor «(k) = J^ 12 (k)z pointing out of the 2D plane. The 
anomalous velocity leads to a net current if k x ft (k) 
has a nonzero average over the electron distribution. In 
the relaxation-time approximation, k is constant over 
the electron distribution. In the parabolic approxima- 
tion, the ground-state distribution of electrons is circu- 
lar. Time-reversal symmetry implies that ft (k) averages 
to zero in the ground-state electron distribution. 

Even without a calculation, we can obtain some intu- 
ition for how the anomalous velocity will lead to helicity- 
dependent photocurrents if ft(k) is proportional to k x , 
which is the symmetry of the GaAs quantum well dis- 
cussed below. Consider a low-frequency circularly polar- 
ized incident wave whose electrical field E lies in the plane 
of the 2D system. If the relaxation time r is short com- 
pared to the period of the incident wave, so that out <C 1, 
then the instantaneous distribution in /c-space is slightly 
displaced from the parabolic minimum, and this displace- 
ment circles with frequency uo around the origin (Fig. 1). 
At an instant when eE points along -fx, the average of 
fl(k) is along +z, while the electron distribution has k 
along ±y according to the sense of circular polarization. 
This gives a current directed along the x-axis, and the 
same direction of current is obtained after a half-period. 

More formally, the current density j arising from the 
anomalous velocity in (Eq. [T]) is 



cPk 
4tt 2 



k x n(k)] fl (k), 



(4) 



where g(k) is the deviation from the ground state con- 
duction electron distribution. Within the relaxation time 
approximation with isotropic mean- free time r, we have 
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eE-v n (k) 

1 jr — iu> 



+ c.c. 



dj_ 

de 



(5) 



and k = eE/h + c.c. (As in the AHE [4], impurity fields 
are not included in the anomalous velocity.) Here v n (k) 
is the normal velocity Vfce(k), and the partial derivative 
of the Fermi function / with respect to energy e is eval- 
uated at the chemical potential. The optical field E is 
specified in phasor notation: E(t) = (E x St + E y y)e %ujt . 
The Berry flux is taken to be ft(k) = /3k x z for some con- 
stant /3 with units of volume, because a Berry flux linear 
in k is required for the effect to appear at quadratic order, 
and any linear combination of k x and k y can be brought 
to the above by a rotation. For a circular Fermi surface 
and low enough frequency that the semiclassical equation 
is valid, the dc current density ] dc is the real part of 
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Switching to polar coordinates, the k integral becomes 



J k 3 dk d0 (cos 0)(E* cos 6 + E* sin 0) (^pj 
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df 
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Ti-^x^-F, 



(7) 



where ep is the Fermi energy. Inserting this result into 
the complex current density (Eq. [6]) gives 



3dc 



f3ne 3 



2h 2 l/r 2 +w 2 



Me x e; - E y E* x )± 



+l/T(E x E* y +E y E* x )±+\E x \ 2 y 



(8) 



In Eq. [H] we substituted cp = 7rh 2 n/m where n is the 
electron areal density and restored a factor of 2 for spin. 

The first two terms on the right-hand side of Eq. [8] are 
Berry-phase contributions to the circular (CPGE) and 
the linear (LPGE) photogalvanic effects, respectively. 
The LPGE, in which the photocurrent is maximal for 
linearly polarized light, is allowed in all acentric media. 
The CPGE, which is maximal for circular polarization 
and changes sign with helicity of the light wave, is allowed 
in the subset of acentric point groups that are optically 
active, or "gyrotropic." The last term on the right-hand 
side of (Eq. is a more typical form of photovoltaic ef- 
fect, in which the asymmetry of the medium fixes the 
sign of the photocurrent, independent of the polarization 
state of the light. But all these effects differ fundamen- 
tally from the photovoltaic effect in a solar cell, for ex- 
ample, where the photocurrent is just determined by the 
number of absorbed photons; here all terms depend on 
the polarization state of the incident light. 



Doped semiconductors are attractive systems in which 
to search for the above effects because of their relatively 
long relaxation times. However, zincblende semiconduc- 
tors have too much symmetry to allow the above effects 
in bulk. Instead, /3 must be generated by a confinement 
potential such as arises in a quantum well structure. We 
will compute this type of Berry phase by two approaches 
and compare the results to photocurrent experiments. 

The first approach is an explicit microscopic calcula- 
tion on model quantum wells. Consider the lowest con- 
duction subband of a (110) quantum well in GaAs and 
assume that the crystal structure (Fig. 2a) is unmodified 
by creation of the well, which then is simply an addi- 
tional electrical potential. The spacing between (110) 
planes is \f2aj\ = 2.00 A with lattice constant a = 5.65 
A. The maximal remaining point group symmetry is C^v'- 
there are two mirrors, in (110) and (110) planes, and a tt 
rotation around [001] that is the product of the mirrors. 
An asymmetric well potential eliminates the (110) mirror 
and the rotation, reducing the symmetry to C s . Either 
symmetry allows a nonzero Berry phase parameter f3. 

Our computational approach starts from reference 
tight-binding parameters for conduction and valence 
bands in GaAs [18]. It is known that this tight-binding 
model gives a roughly correct band structure, but over- 
estimates the band gap and effective mass. Achieving 
agreement with experimental bandgaps and masses re- 
quires correlation methods [19] that are difficult to im- 
plement on a quantum well of realistic size, although the 
envelope expression below explains how these could be 
applied. Tight-binding results for the Berry phase may 
be more accurate than those for the band gap, as the 
former results purely from geometrical properties of the 
wavefunctions. Spin-orbit coupling is ignored since in the 
conduction band it is weak compared to orbital effects. 

The resulting Berry phase factors for a variety of quan- 
tum well potentials are given in Fig. 2b. We find that f3 

o 3 

is on the order of 1 A and decreases with increasing well 
width. For comparison to experiment, consider the cur- 
rent J = j x L from optical power P illuminating a square 
of side L. The CPGE current for circular polarization is 

3mP 
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(9) 



where n was eliminated in favor of the Drude conduc- 
tance G = ne 2 r/m and a = e 2 / (Aneofic) « 1/137 is the 
fine structure constant. The quantity in parentheses can 
be viewed as an effective voltage that combines with the 
Drude conductance. From the above formula, P=l W 
incident on a 1 mmxl mm area of quantum well with 
Drude conductance 10 -2 mho leads to a peak current of 
1.16 nA at cor = 1, where we have taken (3 = 1 A and 
used the effective mass 0.067 m e of conduction electrons 
in GaAs. Currents of this magnitude are found by the Re- 
gensburg group in experiments on doped semiconductor 
quantum wells and heteroepitaxial interfaces [14] [15] [16] . 
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FIG. 2: Berry phase in (110) quantum wells, (a) Positions of 
atoms (Ga gold, As blue) within (110) planes of GaAs show- 
ing inversion asymmetry. Shadowed atoms are in the next 
plane above or below, (b) Strength of Berry phase factor 
(3 = dQ/dk for (110) quantum wells: square wells of depth 
0.25 eV and 0.5 eV, Gaussian well of depth 0.25 eV, and tri- 
angular well of depth 0.5 eV. The wavefunction width 2^/ (x 2 ) 
(horizontal axis), x = (x + y)/\/2 changes with well potential 
and determines (3 almost independently of well shape. 



Extending the tight-binding results given above to a 
broader range of experiments requires understanding the 
confinement-induced Berry phase in an arbitrary quan- 
tum well. We choose an approach similar to the envelope 
approximation for energy bands in a quantum well [20] . 
A three-dimensional zincblende semiconductor like GaAs 
is acentric but nevertheless has no Berry-phase photocur- 
rent at first order in the applied intensity. Microscopi- 
cally, near the Brillouin zone origin T the remaining sym- 
metries restrict the form of the Berry flux to be in the T\ 
representation; Q is not linear in k near T but instead 

^ = A (k X (k 2 y - kl), ky(k 2 z ~ kl),k Z (kl ~ k 2 y )) . (10) 

The constant A determines Berry-phase effects in the ma- 
terial when the Fermi surface is near T. 

A quantum well can reduce the symmetry and combine 
with A to allow a linear photocurrent. Results on the 
(110) well (Fig. 2) suggest that the detailed form of the 
confined state is not too important, so we take a Gaussian 
wavepacket of conduction band states: 



i>2D{x,y,z) = C 



/• 



e i)t - T u k {x,y,z)dk z . (11) 
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This wavefunction combines the intra-unit-cell depen- 
dence u of the conduction band with a Gaussian envelope 
on longer scales. C is a normalization constant. (For sim- 
plicity, we take coordinates with z perpendicular to the 
well plane and later account for the rotation relative to 
the crystalline unit cell.) The wavefunction is localized 
in the z direction, with (z 2 ) — a 2 for a much larger than 
a unit cell, which is the regime where the envelope ap- 
proximation is applicable. 

The above wavefunction is a 2D Bloch state, and the 
2D Berry flux F 2 ® arises from the 3D Berry flux of the 
original crystal. The envelope result is easily stated: the 
2D Berry flux averages the 3D Berry flux over k z . 



F X y (fee? ky) 



xi kyi k z ) dk z 
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dk z 



(12) 



For the (110) well discussed above, inserting the 3D Berry 
phase strength A to compute the 2D Berry-phase coeffi- 
cient (3 gives a simple function of A, the well width <r, 
and a dimensionless factor go set by well direction: 



= 



go A 

<7 2 



A 

"80*' 



(13) 



which is the curve shown in Fig. 2 with A = — 410A 5 
calculated from our bulk tight-binding wavefunctions. 

The same approach to compute the confinement- 
induced Berry phase can be applied to other materials 
and well geometries. For example, the geometrical factor 
of g = -1/8 in (Eq{l3| is modified in a (lln) well to 



1 



9n 



(2n 2 + 4) 3 / 2 " 



(14) 



The vanishing at n = 1 is expected because a (111) quan- 
tum well has an extra symmetry (a 3- fold rotation axis). 

If the Berry-phase mechanism of photocurrents in 
quantum wells can be confirmed experimentally, Berry- 
phase phenomena will have been observed in all the four 
basic materials classes, magnetic and nonmagnetic met- 
als and insulators. However, just as in the case of the 
AHE in TR-breaking metals, other symmetry-allowed 
mechanisms coexist with the Berry-phase contribution. 
Only after many years of debate has it been established 
that the Berry-phase mechanism for AHE is dominant in 
certain regimes of temperature and disorder. 

We believe the Berry phase contribution to photocur- 
rent in acentric metals can be identified more easily, 
through its signature dependence on the frequency and 
polarization state of the light. For example, Eq. [8] pre- 
dicts that the CPGE and LPGE effects are linked; the 
Berry-phase generates both with equal magnitude and 
they vary with frequency like the imaginary and real 
parts of the Drude conductance, respectively [13]. More- 
over, the Berry phase mechanism is readily distinguished 



from alternate mechanisms involving inter-subband ef- 
fects (either spin or orbital) by its dependence on the 
angle of incidence, of the light. While the Berry-phase 
related currents are driven only by in-plane components 
of the optical field, the inter-subband mechanisms are 
driven by E z as well. Thus the Berry-phase photocur- 
rents vary strictly as cos#, while inter-subband response 
has terms proportional to both cos 6 and sin 6. In this re- 
gard, the early measurements of CPGE in a (113) quan- 
tum well by Ganichev et al. [14 are highly suggestive; 
although both cos# and sin# dependences are allowed 
by symmetry, only the cos 6 variation is observed. 

Finally, we suggest that the envelope approximation 
developed here can be applied to confinement-induced 
Berry-phase phenomena in a wide variey of hetero junc- 
tions and quantum wells. Beyond interest in the effect 
per se, measurements of helicity-dependent photocur- 
rents could be used as a means to obtain the leading 
non-zero Berry phase in many materials, in the same way 
as other transport measurements are used to obtain the 
parameters of conventional Fermi liquid theory. 
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